Tuning the shape of the condensate in spontaneous symmetry breaking 
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We investigate what determines the shape of a particle condensate in situations when it emerges 
as a result of spontaneous breaking of translational symmetry. We consider a model with particles 
hopping between sites of a one-dimensional grid and interacting if they are at the same or at 
neighboring nodes. We predict the envelope of the condensate and the scaling of its width with the 
system size for various interaction potentials and show how to tune the shape from a delta-peak to 
a rectangular or a parabolic-like form. 
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The phenomenon of condensation is widely spread in 
nature and appears in equiHbrium systems as well as in 
systems driven far from equilibrium. Examples of equi- 
librium condensation are Bose-Einstein condensation or 
certain phase transitions in quantum gravity In non- 
equilibrium systems, the condensation manifests itself 
in microscopic (intracellular) and macroscopic (highway) 
traffic Q where it corresponds to jamming, granular flow 
[3| and granular clustering or gelation in networks jH| 
where a single node takes a finite fraction of all links. 

Many of these systems can be modeled as a set of 
particles occupying discrete levels, or "boxes" on a one- 
dimensional grid. The balls-in-boxes model (B-in-B) 
or its non-equilibrium version, the zero-range process 
(ZRP) 0, are well-known examples. In these models, 
the condensation transition happens along with a spon- 
taneous breaking of translational symmetry. Although 
these models are more abstract than realistic ones, they 
serve as paradigms since the stationary state can be de- 
rived analytically as fully factorizing over the sites of 
the grid, or, more generally, over the nodes of an ar- 
bitrary graph The factorization is due to ultra- local 
("zero-range") rules assumed to govern the dynamics of 
particles. This makes the phase structure of the system 
accessible and enables comparison with experiments on 
condensation phenomena, at least on a qualitative level. 

As a natural consequence of the lack of interactions 
between particles at neighboring sites, the condensates 
in the ZRP or B-in-B models always occupy a single site. 
The question then arises how the shape of the condensate 
changes in the presence of interactions between sites that 
tend to fiatten out the condensate's profile but still pre- 
serve translational symmetry and conserve the current. 
In this paper we study this problem in a model similar 
to the one proposed in Ref. M. It is related to a solid- 
on-solid (SOS) model (lol . supplied with dynamical 
rules that drive the system out-of-equilibrium. Here, the 
steady state factorizes over pairs of sites which allows for 
nearest-neighbor interactions while making the system 
analytically solvable. We shall show that although for 



a broad class of hopping rates the condensate becomes 
extended to 14^ ~ VlV sites, where TV is the total num- 
ber of sites, it is possible to obtain any scaling W ~ TV" 
with < a < 1/2. We shall also predict the shape of the 
condensate for some special cases. Possible applications 
to surface science will be mentioned in the conclusions. 

We consider a ring with TV sites and M particles of 
unit mass placed randomly on the sites. Each site i can 
carry an arbitrary number rrn = 0, . . . , M of particles. 
The dynamics is divided into two steps. As a first step, 
a particle may leave a randomly chosen site with prob- 
ability u(TOi|mi_i, TOi+i) which depends on the state of 
neighboring sites. As a second step, the particle chooses 
a target to the right site with probability r or to the left 
one with probability 1 — r. One can show fl3| that when 
the hopping rate has the form 



u(mj|mi_i, m^+i) 



g{mi - l,TO,-i) ff(m, - l,m»+i) 
g(m^,mi_i) g{mi,mi+i) 



(1) 

with g{m,n) being a symmetric non-negative but other- 
wise arbitrary function, the steady state factorizes over 
pairs of sites according to 



N 



P(mi, . . . ,TOAr) = Jj5(mj,mi+i)(5 



M 



(2) 

where the ring geometry implies tun+i = mi. The dif- 
ference as compared to the model of Ref. @| is that there 
r — 1 and no assumption on the symmetry o{g{m, n) was 
made. The parameter r changes only the net current of 
particles. When r = 1/2, the current is zero and the sys- 
tem is in equilibrium. Since the steady state ^ does not 
depend on r, its static properties can be calculated by 
formally treating the system as if it was in equilibrium 
with the probability of a microstate given by Eq. |(2]). 

The criterion for condensation was analyzed in Ref. 
by means of the grand-canonical partition function 



Zn{z) 
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where T„ 



,(m+n)/2 



g{m,n) and z is the fugacity. 



Zn{z) must have a finite range of convergence Zc. The 
critical value of the density p = M/N then follows from 



lim 



z d\nZN{z) _ J2,n^'f>l 



N- 



N 



dz 



,2 



(4) 



where (j> is an eigenvector of to the largest eigenvalue 
Amax, for z = Zc- At the critical point, Zjy = A^^^^ for 
large N. In what follows we assume that g{m,n) has 
been rescaled so that Zc = 1, hence the critical density 
can be obtained by diagonalizing Tmn — gim^n). 

To investigate how the introduction of site-site interac- 
tions influences the properties of condensation, we stick 
to the following choice: 



g{m, n) — K{\7n — n\) p{m)p{n) 



(5) 



When K{x) — 1, g{m,n) factorizes and we recover the 
ZRP. Here we assume that both K{x) and p{m) are some 
positive, decaying functions of x and to, respectively. 
The choice ([5]) is motivated by studies on the SOS model 
in the context of surface roughening fiol . and corre- 
sponds to the energy E = — lniir(|TO — n|) — l/2(lnp(TO) + 
Inp(TO)) of the interface in a 1+1 dimensional surface 
(where the surface refers to the envelope of occupation 
numbers). In Ref. the following choice was proposed: 



K{x) 



-Jx 



p{m) 



(6) 



with parameters J and U generating an effective surface 
stiffness and a pinning potential, respectively. For the 
weights we obtain 



a Jo 



(e-^o _ e-2(J-Jo))(e2(J-.Ai) _ l) ' 



(7) 



with Jo = 17 — ln(e'^ — 1). In the limit J Jo, this agrees 
with the asymptotic result in Ref. For J = U — 
1, Eq. ^ gives pc ~ 0.2397, which we also conflrmed 
numerically. If J < Jq, the critical density is infinite. 

The generic case — For the weight functions ([6]), the 
width W was estimated in Ref. [Sl] by a random-walk 
argument to be proportional to y/N. Here we use a dif- 
ferent approach which allows to calculate not only W 
but also the envelope of the condensate in the case when 
K{x) decays exponentially or faster in x, and p(m) 1 
for TO oo. First of all, one cannot use the partition 
function ([3]), because it is not defined in the condensed 
state. We assume, however, that the system can be split 
into the condensate which extends over W sites occupied 
by M' — M ~ N Pc particles on the average and where 
average occupation numbers (to^) grow with TV, and a 
uniform background with (to^) — pc- Since fiuctuations 
in the background are finite and have no long-range cor- 
relations, the mass M' cannot fiuctuate more than ~ ^/N 
so that we treat it as constant. We can therefore assume 



that the probability of having the condensate extended 
to W sites factorizes: 

P{W) « Z^{W)Zc{W) 

cc exp(-M^lnA,nax + lnZe(14^)), (8) 

where Z\^{W) = is the partition function ([3]) for 

the background at the critical point p = Pc and Zc{W) is 
the partition function of the condensate extended over W 
sites and having exactly M' particles. The average exten- 
sion W can be determined by the maximum of Eq. 
Because p{m) 1 for large to, we can neglect the con- 
tribution from p{m) to Zc{W). We have 



w+i 

Zc{W)^ n K{\mk-mk-i 



\fc=i J 



(9) 

where we assumed that nia = mw+i = 0, because oc- 
cupation numbers at the boundaries are small. For large 
systems, the conservation law given by the delta function 
forces the majority of occupation numbers mu to be much 
larger than zero. We can thus extend the summation to 
negative to^ and change the summation over {mk\ into a 
summation over {dk}, where dk = — irik-i- Defining 
the generating function 

oo oo W 

di — — oo d\.v— — oQ k—1 

X s(^-J2kdk-M'^s(^^d,}j (10) 

with auxiliary fields u = {ui, . . . ,uw}, we can rewrite 
Eq. © as Zc(W) ^ G(W + 1,0). The first delta func- 
tion in Eq. lfTO|) implements the mass conservation in the 
condensate, the second delta refiects fixed-boundary con- 
ditions. Using integral representations of both deltas, the 
generating function l|10p can be evaluated in the saddle 
point as G{W,u) ~ e^^^'"'"\ where 

w 

F{z,v,u) = -M'z + Y^^^i'^k + V - kz). (11) 

k=l 

The variables z — z{u), v ~ v{u) are determined from the 
saddle-point equation dzF{z,v,u) — dyF{z,v,u) = 0, 
and the function K{x) is defined by 



K{x)^ J2 Ki\d\)e^'. 



(12) 



We assume that this sum exists, as it does for the case 
From the definition of G{W, u), (dk) = -£^F[z, v, 'ii)|ff=o- 
The symmetry of the averaged peak implies that (dk) = 
— (dw-k) and hence z — {2/W)v. We obtain 



(dk) = k'ix)/Kix) 



x=v{l~2k/W) 



(13) 
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In addition, from the conservation law J2k ^ i'^k) 
we have for large systems 



1 
2^ 



" xK'jx) 

k{x) 



Ax 



-M' 



(14) 



where we defined the reduced extension w = W/^jM' . 
Now, since hiZc{W) = F{2v/{wy/M'),v,d), we can 
rewrite Eq. ([HI) as a function of v: 

^ = -wlnAmax + - / \nK{x)dx, (15) 



having in mind that w is not an independent variable 
but is bound to v through Eq. (fT4l) . Taking the deriva- 
tive with respect to v we find after some calculations that 
P{W) has a maximum for vq = K~^{Xmax)- The corre- 
sponding width W = Wo V M' of the condensate can be 
determined from Eq. lfT4|) by calculating wq numerically. 
Since wq is independent of the system size, we have thus 
shown that W '-^ VM' '-^ ^fN . A detailed discussion will 
be presented elsewhere [l^ . 

To calculate the envelope of the condensate, it is con- 
venient to define h{t) = {nin) / \JM' in a rescaled variable 
t — — — 1, which removes the dependence on the sys- 



tooVM' 

tem size. In Fig. [T] we show that the Monte Carlo (MC) 
data points, obtained for different sizes, collapse onto a 
single curve that looks like a parabola. To calculate it 
analytically, we note that (m„) = X]fe=i ('^fe)- Changing 
the summation into integration, we obtain 



2vo K{vQt) 



For the special case of Eq. ([6]) , one obtains 

h{t) 



Wq , /cosh J — cosh Wo * 

in 

2vo \ cosh J — cosh vq 



(16) 



(17) 



with Wo, Vq being functions of the parameters J, U. The 
formula for vq is particularly simple, vq = J — Jq. The 
comparison between Eq. lfT7|) for J = J7 = 1, which gives 
vo = 0.5413 and wq = 2.2005 (from Eq. and MC 

simulations in Fig. [1] shows a very good agreement. 

Extension different from ^/N — So far we consid- 
ered the situation which, in the language of the SOS 
model, corresponds to interactions between neighboring 
sites which decay rapidly with the difference of heights, 
and the pinning potential Inp(m) localized near zero. 
In this case, the extension scales generically as ~ \/N. 
We will show now that the exponent a in the extension 
W ~ can be tuned to a ^ 1/2 for the following choice 
of the weight functions: K{x) ^ e~^^ and p(m) ~ e"™"' 
with /3,7 > 0. Such weights translate to the SOS en- 
ergy E = \mi — mi-i\^ + m] . The exponent 7=1 



corresponds to a constant, e.g., a gravitational field [ll| 
acting on the envelope, while 7 < 1 and 7 > 1 corre- 
sponds to some attractive force, decreasing or increasing 




Figure 1: Comparison between the envelope of the condensate 
from Eq. I|17p and rescaled MC data for the weights ^ with 
J = U = 1 {or N = 1000, p = M/N = 1 and 3 (circles, 
squares) and A*' = 4000, p = 1, 3 (diamonds, triangles). 




Figure 2: (color online) Phase diagram for K{x) ~ e ^ and 
p{m) ~ e"™^. The exponent Oroct = (/3 - 7)/(/3 - 7 + 1) 
for rectangular and Osmooth — {l3 — 7)/(2/3 — 7) for smooth 
condensates. The dotted lines show a = 0.05, 0.1, 0.45 
(from left to right). 



with m, respectively. Although the series (fT2l) is con- 
vergent, Imim-yoo Pim) — > 0, so that we cannot use the 
previous method. However, we can proceed as follows. 
From the condensation criterion we obtain that for 7 > 1 
there is no condensation. The condensate emerges above 
some pc only for 7 < 1, the value of (3 does not matter. 
If the condensate was localized at a single site, its sta- 
tistical weight would be Pi ^ Nc^-^p{N)K^{N). Here 
c stands for the contribution of background sites. For a 
two-site condensate, both sites will be almost equally oc- 
cupied because every difference e is suppressed as if(e), 
and we have P2 - Nc'^-^K^ {f)p^ (f ). We thus ob- 
tain lnPi/P2 ~ —cpN'^ + CyN'' with some constants 
C(3,Cj > for 7 < 1. This means that for 7 > /3, the 
ratio P1/P2 ^ 00 as — > 00 and the condensate is lo- 
calized. However, for -f < (3, the contribution from the 
two-site condensate is larger. Similarly, one can show 
that the contribution P„ for an n-site condensate grows 
with n, so the condensate must be extended. Denoting 
by P{W) the probability for an extension W, we can 
estimate \nP{W) as follows, dropping inessential con- 
stants. First, there is a negative — W term that accounts 
for the W excluded background sites (cf. Eq. An 
analogous term comes from fiuctuations (we assume that 
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the entropy is additive). Second, every site in the con- 
densate contributes as p{N/W) (because M' ~ A''), so 
there is a —W{N/W)'^ term. The third term comes from 
the weight factor K{x) and depends on the condensate 
shape. As originally suggested by numerical experiments, 
we only have to distinguish two types (see below). For 
shapes with a smooth envelope h{t), the contribution is 

r^W j \nK{2Nh\t)/W^)dt - -W (N/W^)'^ , (18) 

while for a rectangular shape, where h{t) has a disconti- 
nuity at the borders, the contribution is 

2lii K (N/W) ^ ~ {N/W f . (19) 

Taking all terms together, the extension W follows by 
searching the maximum of the larger of the two values 
P{W): 

In Psmooth ^-W- W{N/Wy - W {N/Wy , (20) 

InPrectangular ^ -W - WiN/Wy - [N/Wf . (21) 

Assuming W ~ N", this gives arcct = (/? — 7)/(/3 — 7-t- 1) 
for the rectangular, and asmooth = (/3— 7)/(2/3— 7) for the 
smooth condensate, respectively, a non-trivial exponent 
that reminds of finite-size scaling at second-order phase 
transitions. Mixed shapes, partially smooth but with 
some discontinuities, would give a contribution which is 
always smaller than the larger value of Eqs. (|20|) and 
(|2T|l . Thus we can have only these two types of extended 
shapes, and a localized one, depending on the values of 
P and 7. In Fig. [21 we present the full phase diagram. 

A similar reasoning can be applied to other weight 
functions. For instance, for K{x) cx x~'^ (long-range in- 
teractions) and p{m) = e^"^™" , the corresponding proba- 
bilities are 

InPsmooth r^-W^ ,yWln{N/W^), (22) 

InPrectangular ^ -W - l^HN/W), (23) 

and the latter probability is always larger, provided that 
W « const, so the shape turns out to be always rectan- 
gular. Similarly to finite-size scaHng at first-order tran- 
sitions, the height scales with N, while the width stays 
approximately constant, corresponding to a correlation 
length for first-order transitions that stays finite. For 
another choice, K{x) oc x~'' and p{m) cx m~^, one can 
show that the condensate emerges only if 6 > and 
ly > 1 — 6/2, as follows from Eq. (0]), applied to the eigen- 
vector (pm ~ m^''/^^''. If the condensate exists, it is 
always localized, because the ratio P1/P2, calculated as 
before, behaves as ^ N'' and grows to infinity for 6 > 0. 
Even if one performs the sum over the difference e in oc- 
cupation numbers. Pi is still larger than P2 by a factor 
of at least N^^'^. Details are postponed to Ref. [l3|. In- 
terestingly, from experimental observations of the shape 



and the scaling of the width with the system size, one 
can trace back the class of hopping interactions that are 
compatible with the observations. 

To summarize, we have shown that a condensate which 
emerges in a PFSS as a result of spontaneous symmetry 
breaking, can be either extended or localized, depending 
on the competition between local and ultralocal inter- 
actions, and that its extension and the envelope can be 
calculated analytically. Let us close the paper with a re- 
mark on possible appHcations. When atoms condense on 
a crystal surface, they can migrate and build extended 
islands. As experiments on the deposition of clusters [3| 
or fabrication of quantum dots [141 show, the islands can 
be extended in the direction perpendicular to the sur- 
face. It would be interesting to check if the PFSS model 
could be used to predict shapes and typical sizes of is- 
lands of atoms obtained in this way. Also problems like 
mass transport on arbitrary networks could be addressed 
within the PFSS formalism. 

B.W. and W.J. thank the EC-RTN Network "EN- 
RAGE" under grant No. MRTN-CT-2004-005616for sup- 
port. B.W. thanks the International Center for Transdis- 
ciplinary Studies (ICTS) at Jacobs University for support 
of some visits during this collaboration. 



[1] P. Bialas, Z. Burda, B. Petersson, and J. Tabaczek, Nucl. 
Phys. B 495, 463 (1997); Bas V. de Bakker, Phys. Lett. 
B 389, 238 (1996). 

[2] D. Chowdhury, L. Santen, and A. Schadschneider, Phys. 
Rep. 329, 199 (2006). 

[3] C.-H. Liu et aL, Science 269, 513 (1995); S.N. Copper- 
smith et aL, Phys. Rev. E 53, 4673 (1996). 

[4] S. N. Majumdar, S. Krishnamurthy, and M. Barma, 
Phys. Rev. Lett. 81, 3691 (1998); J. Stat. Phys. 99, 1 
(2000). 

[5] P. L. Krapivsky, S. Redner, and F. Leyvraz, Phys. Rev. 

Lett. 85, 4629 (2000). 
[6] P. Bialas, Z. Burda, and D. A. Johnston, NucL Phys. B 

493, 505 (1997). 
[7] M. R. Evans, Braz. J. Phys. 30, 42 (2000). 
[8] M. R. Evans, S. N. Majumdar, and R. K. Zia, J. Phys. 

A: Math. Gen. 39, 4859 (2006). 
[9] M. R. Evans, T. Hanney, and S. N. Majumdar, Phys. 

Rev. Lett. 97, 010602 (2006). 
[10] S. T. Chui and J. D. Weeks, Phys. Rev. B 23, 2438 

(1981); T. W. Burkhardt, J. Phys. A: Math. Gen. 14, 

L63 (1981). 

[11] J. M. J. van Leeuwen and H. J. Hilhorst, Physica A 107, 
319 (1981). 

[12] B. Waclaw, J. Sopik, W. Janke, and H. Meyer-Ortmanns, 
in preparation. 

[13] C. Brechignac et al., Phys. Rev. Lett. 88, 196103 (2002); 

B. Yoon et aL, Surf. Science 443, 76 (1999). 
[14] R. Notzel, Semincond. Sci. TechnoL 11, 1365 (1996). 



